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I. INTRODUCTION 



o : 
o . 
o ■ 

' One of the outstanding issues in gravitational wave research is calculating the wave output from the last stages of 
1 inspiral of binary black holes (BBHs). This intermediate binary black hole (IBBH) problem has been discussed by 
D Brady, Creighton, and Thorne The purpose of this paper is to provide an approximate four-dimensional BBH 

metric from which initial data can be extracted and evolved numerically into and through the IBBH region, 
^f-^ ■ The approach I take is based on the work of Manasse j|] and D'Eath [^|J4|. I consider two widely-separated non- 
\ rotating black holes in a circular orbit. The black holes' mass ratio is not restricted — thev can have comparable masses. 
C\l ■ However, the masses are assumed to be much smaller than the distance between therry. As a result spacetime can be 
^ divided into four regions, each with its own approximation scheme to solve Einstein's equations. There is a strong- 
. gravity region near each of the black holes which is described by the Schwarzschild solution plus a perturbation due 
' to the companion's tidal field. This perturbation is constrained to satisfy the linearized Einstein equations (LEEs) 
, about the Schwarzschild metric. The companion black hole's electric-type and magnetic-type tidal fields are both 
' taken into account in calculating the perturbation. 
^\ . Outside the strong-gravity regions but within the near zone, the metric can be approximated by a post-Newtonian 
On ■ expansion. Further out is the radiation zone which contains outgoing gravitational waves and can be described by a 
[ post-Minkowski expansion of the metric. 

There are overlap zones in this spacetime where the regions described above intersect in pairs. In the overlap zones, 
^ two different approximation schemes — one from each of the two intersecting regions — are both valid. The perturbative 
(2J[), expansions produced by the two approximation schemes are matched in the overlap zones using the framework of 
1^ • matched asymptotic expansions. The post-Newtonian near-zone metric — taken from Q — and the radiation-zone 
. ^ ' metric — taken from Q — already match in their overlap region. In this paper, the post-Newtonian near-zone metric is 
matched to a perturbed Schwarzschild metric in the matching or buffer zone surrounding each black hole. This yields 
information on the asymptotic behavior of the Schwarzschild perturbation at large distances from the horizon, and 
on the coordinate transformation between the two buffer-zone coordinate systems. 

The Schwarzschild perturbation and coordinate transformation are not uniquely determined. However, a different 
choice of transformation — and hence different form of Schwarzschild perturbation — should still represent the same 
physical situation. In other words, different perturbations that match to the post-Newtonian near-zone metric are 
expected to be related via gauge transformations. For the purposes of this paper, it is sufficient to find one trans- 
formation and one Schwarzschild perturbation associated with each black hole that result in a match between the 
post-Newtonian near-zone metric and the distorted-black-hole metrics. 

An approximate spacetime metric is put together by joining the regional metrics at some specific 3-surfaces in the 
matching zones. The final 4-metric is written in a single coordinate system valid up to (but not inside) the apparent 
horizons of the black holes. This metric is useful not only as a source of initial data for numerical evolution, but also 
as a check on the early stages of such an evolution. 

It has been suggested that numerical simulation of BBHs should be performed in corotating coordinates |^ . For 
this reason the metric in final form is given in corotating coordinates. The BBH spacetime can be sliced and spatial 



^Throughout this paper I use geometrized units in which G — c — 1. 
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coordinates chosen in any convenient way when extracting initial data from the metric. (Asymptotically inertial 
coordinates can be used, for example.) 

Initial data generated by the method presented in this paper have the advantage of being connected to the early 
inspiral phase of the BBHs. Detailed gravitational waveforms from this early inspiral phase have already been 
calculated using post-Newtonian expansions. These waveforms will be easily linked to the waveforms obtained by 
evolving initial data extracted from the metric presented here. 

Initial data from this metric have the additional advantages of not containing spurious gravitational waves and of 
reliably describing the physical situation of coalescing BBHs. The accuracy of this description can be improved by 
taking the calculation in this paper to higher orders. 

In Sec. |l|, the near-zone and radiation-zone metrics are written down. In Sec. HE the first black hole's tidal 



deformation is calculated. In Sec. IV, the buffer-zone coordinate transformations are determined, and the distorted- 



black-hole metrics are written in corotating post-Newtonian coordinates. The full spacetime metric is summarized in 
Sec. 0. 

II. NEAR-ZONE AND RADIATION-ZONE METRICS 

Blanchet and collaborators ( and references therein) and Will and Wiseman have calculated in detail the near- 
zone and radiation-zone gravitational fields of compact binary systems. The approach taken by Will and Wiseman 
is particularly useful here because they use a single coordinate system — harmonic coordinates — to cover both the 
near zone and the radiation zone. As a result, expressions for the radiation-zone metric components taken from |^ 
automatically match (to some finite order) the harmonic-coordinate, post-Newtonian, near-zone metric components 
calculated in Q. For this reason I work initially in harmonic coordinates {t' ,x' ,y' , z') with the origin of the spatial 
coordinates placed at the binary system's center of mass. I use only the first post-Newtonian (IPN) metric, not the 
full 2.5PN metric given in j5||^ Consistently with this, I put the black holes on Newtonian trajectories: they are 
taken to be in circular orbits with Keplerian orbital angular velocities. Moreover, I use the post-Newtonian metric 
for point-like particles; in the near zone, I ignore the black holes' internal structure. The near-zone gravitational 
effects of the black holes' multipole moments can in principle be computed by matching out to the near zone the 
tidally-distortcd Schwarzschild metrics obtained in this paper. However, these effects are too small to be included in 



this paper; this is discussed further in Sec. iVB 



A. Binary-system parameters 

Label the black holes BHl and BH2, and let mi and TO2 be their respective masses. Define 

m = nil + "12, ora = nii — m2, /i = . (2-1) 

TO 

Denote the harmonic-coordinate trajectories of the black holes by x\{t') for A ~ 1,2 and j = 1, 2, 3. In other words, 
Xj^{t') are the spatial coordinates at time t' of the center of attraction of the gravitational field of black hole A. 

In this section, boldface letters are used to denote spatial coordinates. For example xa = ix\,x\,x\) = 
(xA, i/A, za)- The notation a • b is used for the quantity Sjka^h^, and |a| is by definition (a • a)^/^. 

Denote the black holes' separation |xi — X2I by h. The circular, Newtonian trajectories of the black holes are 

xi(0 = — b(0, X2(0 = -— b(t') (2.2) 

TO TO 

where 

b(t') = xi(t') - X2(i') = 6(coswt',sintjt',0) (2.3) 

and 



^Higher order versions of this calculation will presumably use higher order post-Newtonian metrics. 
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is the Keplerian orbital angular velocity. Define 



for A = 1,2. By assumption, e ^ 1. 



J, r=(x'2 + y'2 + z'2)i/2^ 



= |x - x^l 



dt' ' 



rA 



V = vi — V2 = e(— sinwt', coswt', 0), 



(2.5) 



B. Demarcation of four regions in the BBH spacetime 



Let us first fix precisely four regions in this BBH spacetime; each of these regions will receive a metric calculated 
as an approximate solution to the Einstein equations. With such a partition of spacetime in mind, define the inner 
limits r™ — yfrn^h and r™ = \fm^. These arc just convenient choices for the inner limits. The important property 
r'l" has is that both r'i"/6 and mi/r'^" — > as mi/6 0. Similarly r|'7^ ^ and m2/r^ ^ as 7712/6 0. Also 
define the outer limit r°"* = Ac/27r — &/2e where Ac = tt/lo is the characteristic wavelength of gravitational radiation 
emitted by the binary system. 

Divide spacetime into four regions that are bounded by the black holes' apparent horizons and the surfaces ri = r™, 



r2 



and r 



^out . 



(i) the region ri < r™ (but outside the apparent horizon of BHl), labeled region I; (ii) the 
region r2 < r™ (but outside the apparent horizon of BH2), labeled region II; (iii) the subset of the near zone specified 
by Ti > r™, r2 > r™, and r < r°"*, labeled region III; and (iv) the region r > r°"*, labeled region IV. The near zone 
contains region III and overlaps with regions I and II; the radiation zone corresponds to region IV. The buffer zone 
around black hole A contains ta — and satisfies itla <C r^i <C 6. These regions of spacetime are illustrated in 
Fig. I 




FIG. 1. Schematic illustration of the various regions in the BBH spacetime. Regions I, II, III, and IV are demarcated by 
solid lines; the buffer zones are bounded by dashed lines. 
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C. Near-zone metric in harmonic coordinates 



In the near zone, the IPN harmonic-coordinate metric with two point-hke particles representing the black holes is 

1 

go'O' =-1 + - + ^-2^ + ^ +^ [Avl - ni • vi 2 + -1 [Avl - n2 • 

ri r2 \ ri r2 / ri '- ' r2 ^ 

mim2 / I 1 \ mim2. , . 

-b • (ni - n2), 



Tl 7-2 



, , mi „ rrii ■ 
' r-i r2 



53 



, 2toi 2m2\ 

5,,,, =,5,^ ( 1 + — i + — ij . (2.6) 



7'1 ^2 



This metric presumably differs in the near zone by a small amount from an exact solution to the Einstein equations 
representing BBHs. I take the neglected terms in the 2.5PN metric P| to be an estimate of the errors in the IPN 



metric (2.6). 

The largest neglected terms in (7o'0' are of the form rr? jlP'TA-, m? /hr\, m\/r\, m^/b^, and errv'r^ /b^. (The last 
term represents a radiation reaction potential.) Let us compute the orders of magnitude of these terms at various 
places in region III. If > r™ ^ be (here and henceforth "~" means "is of the order of" and a > & means a> b and 
a ^ b) for A = 1 or 2, then the error in go'o' (denoted SgQ'o') is of 0{e^) and comes from neglecting a term of the form 
m\lr\. If both ri ^ b and r2 ~ 6, then 5go'0' ~ e^- Finally, if r ^ r""* ~ 6/e (so that ~ 6/e for ^ = 1 and 2), then 
the error (5go'0' ^ arises from neglecting the radiation reaction potential. Note that it is reasonable to consider the 
"absolute" errors Sg^i^' in the metric components since the coordinate system being used is asymptotically inertial 
and the errors are only calculated in regions of weak gravity where deviations from a flat metric are small. 

A similar analysis for go'i' yields (Jgo'i' ^ if fA ^ for A = 1 or 2, 8gon' ^ if both n ~ 5 and r2 ^ b, and 
Sgo'i' ~ if 7" ^ ^out^ Lastly, Sgi'j' ~ if rA ^ t"™ for A = 1 or 2 (this comes from neglecting a term of the form 
TO^ /r^ in gi'jr ) , Sgi'j' ~ if both ri ~ 6 and r2 ~ b, and Sg^'j' ~ if r ^ r""* . 



D. Near-zone metric in corotating coordinates 



The metric (2^) is transformed to corotating coordinates (t, x, y, z) defined by 

t' — t, x' = X cos ujt — ysincjt, 
y' — X sin Lot + y cos Lot ^ z' = z. 

In terms of the new coordinates, 

r^{x^ + y^ + z^Y'\ 

Putting the expressions (2^) - (2^) in Eq. ( ^.(j ) and transforming to corotating coordinates gives 
,™ ,™ , ™ ™ , 2 

ds^ = dt^' 



(2.7) 
(2.8) 



2mi 2m2 
-1 + - + --2 



ri 



r2 



mi m2 

ri r2 



3/i / m2 mi 
b \ri r2 



1 



1 



n r2 



1 



2mi 2m2 



8/ie ( — - — ) dtdy + ( 1 



2mi 



f2 
2to2 



(.2+y2) 



2lo I 



mi 
~3" 



1 1 

7-1 7-2 



2mi 2m2 \ 
^2 / 



dt{xdy — ydx) 



(2.9) 



(dx^ + dy^ + dz^), 

fi r2 y 

where, in terms of the new coordinates, the quantities ta are 

ri — [{x — 'm2b/rnf' + + z^] ^^'^ . 

r2 = [{x + mib/m)^ + + ■ 

This is the final form of the metric in region III. (Note, however, that this metric is valid throughout the near zone, 
which includes the buffer zones around the black holes). It remains to specify the metric in regions I, II, and IV. I 
postpone until Sec. IIP discussion of the errors Sgfn, in the metric components in the new, rotating coordinate system 
{t,x,y,z). 



(2.10) 
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E. Radiation-zone metric in harmonic coordinates 



The radiation-zone metric can be extracted from In that paper, a potential /i'' ^ is defined by 



(2.11) 



where rf" " =diag(-l, 1,1,1), Q^i'v' is the spacetime metric, and g' = det((/^tV')- Equation (5.5) of gives ^ in the 
radiation zone in harmonic coordinates {t' , x' ,y' , z') for a system of several bodies. After correcting a typo in that 
equation]^, I specialize to a system of two bodies of masses mi and rn2 in a circular orbit specified by Eqs. (2.2)-( p75| ). 
This yicldsR 



h° " {f, x', y', z') = + + 2 



/i°'''(t',x',2;',z') 



r r 
1 



r 



Q'J [u') - e'^^ J\u') 



r 

-I- 

3 1 r 



,ijk 



Q'i''[u')-2e'^^f^{u') 



where 



and 



h''^'{t',x\y',z') = ^n"n'^" + -Q'^u') -U- \q^^' {u') ~ Ae'^^^'^' j'^^Hu') 

r 6 \ r I 



rh — m(l — ^/26), u — t' — r, n' — x'/t" 



(2.12) 



(2.13) 



Qij ^ ^ mAX^x^ = fibV, Q""^^ ^ ^ mAXAx'^x'X = -fi{Sm/m)bVb'', 

A=l A=l 
2 2 
7^ \ ^ _ ilm I m ilmil m jij \ ^ _ ilm I mi /' ilmil mij 

J = m-A^ ^A^A ~ 1^^ V , J ■' ^ y. "^Af ^A^A^A ~ ~IJ.[om/m)e o v w . 



A=l 



A=l 



(2.14) 



Putting the expressions (2.14) in Eq. (2.12) and using Eqs. (|2.2|)-(2.5) gives 



'™ ' % + ^ (2(n' • v)^ - |^(n' . b)^ + ^(n' . b)(n' • v) + -L [3(n' • b)^ - 5^] 



^^|Zj(n'.b)^(„'.v)-2(„'.v)3 + l(n'.b) 



^(n'.b)^-12(n'.v)^-^ 



+ -|(n' • v) [6^ _ 5(n' • b)^] + ^(n' • b) [ih^ - 5(n' • b)^] |, 



(n'-v) + -(n'-b) 



- S(n' . b)64 - ^ — (-S(n' • b) [3(n' • hy + 4(n' • v)6^] 
J r m O'' 



+ 2(n'-v)"w* + - 

r 



6(n' • b)(n' • - ^(n' • b)^^^ + "^b' 
¥ b 



4 [3(n' • b)^ - b'] 



1 



b-^ 



2^ (5m f 6m , 

+ — <^ -^{n' ■hy'V> 

r m \ 



(n'-v) + -(n'-b) 

r 



where v and b are evaluated at the retarded time u' = t' — r. 



(2.15) 



^The term Am/r' in the expression for /i"*^ should instead be Arh/r' . 
*Note that I have replaced r' in Eq. (5.5) of M by r. 
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The metric g^,'!/' can be gotten from Eq. (2.15) as follows: from Eq. (2.11) we have 



(2.16) 



Take the determinant of both sides of Eq. ( 2.16 ) ; this yields g' — det (r/^ — h^^ ). So g' can be calculated once h'^ ^ 
is known, and then ^ can be gotten from Eq. ( 2.16| ). Inverting the matrix g^ gives the spacetime metric g^'yi . 

When performing these calculations, I keep all terms of the form to3-p/2^-3(i-p/2)^-p integer p > 0. I also 
keep — at each order in r — all terms that are of lower order in m/b than this, and throw away terms of higher order in 
m/b. This means in particular that no terms of 0{r~^) are kept. This scheme of organizing terms is consistent with 
the ordering of terms in Eq. (5.5) of 

The result of these calculations is the following radiation-zone metric in harmonic coordinates: 



2™ 2to^ /i 



fi 6m 



{(n'-v) 



r 
7m 



2(n' . v)2 - ^(n' • b)2 + V' • b)(n' • v) + 1 [3(n' • b)^ - b^] 



-(n'.b)^-2(n'.v)^-:^ 



+ -(n'.b) 



3m , , , , 9 , , , , TO 
-^K-b)^-6(n'-vf-- 



^(n' • v) [b^ - 5(n' • b)2] + l(n' • b) [3b^ - 5(n' • b)^] 



4m 

go'i' = 

r 



1 



(n'-v) + -(n'-b 



^(„' . b)(n' . V) + 1 [3(n' . hf - b'] y + |-^(n' • b)(n' • v) + ^ 



g,,y = 1 



2to m^ M r « / / N 9 2to ,,,,9.6,,,,,, , 1 



u 5m ( . , , 
r TO [_ 



7m 



2(n' • v)2 - -^(n' • b)2 + ^(n' • b)(n' • v) + [3(n' • b)^ - b'] 



^3 (n'.b)^-2(n'.v)^ + ^ 



^(n'.b) [^(n'.b)2-2(n'.v)2 



A(n' . V) [b^ - 5(n' . b)^] + l(n' • b) [3b^ - 5(n' • b)^] }^ + ^n'^n'^ + ^ (v^v^ - ^bV 



2/i Sm ( 6to , , , , 
r m b -^ 



(n'-v) + -(n'-b) 

r 



/ m 



(2.17) 



where v and b are evaluated at the retarded time u' = t' — r. 

The errors dgfj_i^i in these metric components in region IV can be estimated by computing the orders of magnitude 
of neglected terms, which are of the form to^~p/^6~^'^^~^/^V~''(to/&)"/^ for integers p > and n > 0. These terms 
include m^/b'^r, eTO^/r^, m'^b/r^, etc. This gives 5g^i^i ^ t' for r > r°"* ~ 6/e, 55^',^' ~ for r ~ bjt^ ^ j.o«t^ 
Sgf,'^' < for r > b/e'^. 



F. Radiation-zone metric in corotating coordinates 



Substituting the expressions (2.2)-(2.5) into the metric (2.17) and transforming this metric to corotating coordinates 
{t,x,y,z) [defined in Eq. (2.7)] gives 



2m / /i \ 2m^ 2e , ^ 

— 1 H 1 — h A H — — (xcoscjr — y sm ujrjB — 2e(a:smwr + y cos ujrjD 

r V 26/ r-^ b 



+ uP'lx^ + y'^)E H — ^(a;costi'r — ysina;r)^A^ + e^(a;sinijjr + ycoswr)^^ 



12/Lte^ 5m 



(xcoswr — y sina'r)'^(a;sina;r + y coswr) 



r^fe^ m 

2dt \ €[sui{ijjr)dx + cos(ajr)d?/]_B + b[cos{ijjr)dx — s\n{ujr)dy]D + uj{xdy — ydx)E 
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-{x COS Lor — y siiiLur)[siTi(Lur)dx + cos{Lur)dy]N 



where 



+ 



eb{x sinujr + y cos LL!r)[cos{ujr)dx — sm{u!r)dy]S 



r'^b m 



(x cos ujr — ysinujr) [(xcos2a;r — y sm2u}r)dx — (a;sin2a;r + y cos2ujr)dy] 



+ E{dx^ + dy^ + dz^) H -rixdx + ydy + zdzf' + echini ujr)dx + cos{ujr)dyf N 

+ b'^[cos{ijjr)dx — sin(ti)r)dy]^5 
Xlfji^ Sui 

H 2 (a;coswr — y siiiujr)[cos{LL!r)dx — sm{ujr)dy][siii{ujr)dx + cos{ujr)dy] 



(2.18) 



A 



D 



2e^ [(a;sincL;r + ycosLur)'^ — {x cos cor — j/sinwr)^] 



H [xsmLjr + y cosujr)(x cosuir — ysmujr) H — j |3(a;cosa;r — ysmujr) — r 



^ Sm ( b^ 



-(xcoscor — ysinujr) [3r^ — 5{x cos ujr — ysmcjrY 



3bh 



H — (x sin Lor -\- y cos ujr) I r — 5(x cos Lur — y sin urY 



2m 



(a;cosa;r — y sinojr) [3(a;cosCiJr — ysm 



ojr) 



6{x sinwr + y coswr)^ 



+ e^(a;sinajr + ycoscur) [7(a;cosu;r — j/sinwr)^ — 2(a;sina;r + j/coswr)^ — 



5 = -^ 



e(x sin ur + y cos ujr) H — (x cos Lor — y sin ldt) 



2/i 5m for/ \2 i \2^ 

<e \2ixsniLor + y cos Lor) — iixcosLor — ysmur) 

r'^ TO 

^ — (xsinwr + y cos wr) (a; cos — y sinccr) H — - [3(a:coswr — y sincjr)^ ^ j^' 



2^16^ 



- [x cos Lor — y sin Lur) 



1 



4e 



5m ( 
r m \ b 



(x cos Lor — y sin ) (x sin tor + y cos cur) 



E ■ 

N 

S 



H — [3(xcoswr — ysintjr)^ — r^] 



2to / ii\ m" 



2^e^ 5m 



- {x cos u)r — y sin ojr) + e (x sin wr + y cos wr ) 



r 

24xe^ 
rb 



5m 

TO 



- [X sm ujr + y cos ujr) H — ^ (x cos tor — y sm j 



2 1 5m 
6 r TO 



- (x sm uir + y cos tor) H — (x cos lot — y sm j 



(2.19) 



and e = (to/6)^/^, w = [m/b^Y^''^ . This is the final form of the metric in region IV. 

It is now necessary to evaluate the errors 5g^i, in the metric components in the corotating coordinate system 
(i, X, y, z). Since this coordinate system is not asymptotically inertial, it no longer makes sense to compute absolute 
errors. The rotation of the coordinates introduces terms of Oiujr) and 0[(tjr)^] in goo and terms of 0{ujr) in g^a. 
For this reason, I define the "normalized" errors J^qo = 5goo/ {'-orY , Sg^^ — 5gQi/u)r, and 5g^^ = 5gij in region IV. It 
follows that Jg^^ ^ e' for r > r°"* b/e, Jg^^ for r ^ b/e^ > r""*, and Jg^^ < for r > b/e^. 

In region III, ojr is less than 1 (a nd in the buffer zones, wr <C 1). So rotation of the coordinates is not important in 
analyzing errors in the met ric (2.9) in region III. I continue to use absolute errors in that region. The errors Sg^j^i, in 
the metric components (|2.9|) in corotating coordinates in region III are the same as the errors in harmonic coordinates 
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(see Sec. II C ): (i) Sgoo e"^, Sgoi ~ e'^, and Sgij ~ if va ?J for ^ = 1 or 2; (ii) Sgoo ^ e^, 6goi ~ e^, and Sgi 
if both ri ~ 6 and r2 ~ b; and (iii) ^g^i^ ~ if r < 7.°"*. 

Since the analysis by Will and Wiseman ||^ of con ipac t binary systems uses a single coordinate chart to cover both 



the near and radiation zone s, the near-zone metric (2_^) automatically matches (to some finite order; see below) the 
radiation-zone metric (2.18) at r = r°"*. The match is not perfect because I have truncated the relevant perturbative 
expansions at finite order. As a result, there are discontinuities in the metric components at r ~ 7-0"*. The orders 
of magnitude of these discontinu ities can be estimated as fo llows: first expand in powers of b/r for r > b and 



substitute this expansion in Eq. (|2.9|); then expand Eq. (2.18|) in powers of Lor for r < r°"*; finally, compare the two 



The result is that the discontinuities in 5^^, denoted [g^v]^ are [gf^u] 



III. TIDAL DEFORMATION OF THE FIRST BLACK HOLE 



The metric (2.9) is valid not only in region III but also in the buffer zones around the black holes. The next step is 
to match this metric to a tidally-distorted black-hole metric in the buffer zone around BHl. There are two coordinate 
systems which ove rlap in the buffer zone. The first is the corotating post-Newtonian coordinate system {t,x,y,z) 



defined in Eq. (2/7). The second — to be called the internal coordinate system — covers the strong-gravity region near 
the first black hole and is valid from the black hole's apparent horizon up into (and through) the buffer zone. The 
internal coordinates are chosen to be isotropic coordinates (T, X, Y, Z) in which the unperturbed Schwarzschild metric 



where 



The region these coordinates cover will be called the internal region; it contains region I in particular. 

The first step in matching the near-zone metric (|2.9| ) to the internal metric [which is the Schwarzschild metric (3.1) 
plus tidal perturbations] is to write the metric ( |2.9|) in internal coordinates. Then the near-zone metric and the 
internal metric are both expanded in positive powers of mi/R and R/b in the buffer zone. Finally, corresponding 
terms in the two asymptotic expansions are equated. The near-zone metric determines in this way the asymptotic 
form of the tidal pert urba tions on BHl. These perturbations arc further constrained to solve the LEEs about the 
Schwarzschild metric ( |3.l[ ) and to be finite at the horizon R — toi/2. 

The asymptotic form of the Schwarzschild perturbations in internal coordinates can be determined — independently 
of the matching procedure described above — by calculating the electric- and magnetic-type tidal fields of BH2 in the 
buffer zone surrounding BHl. Once this asymptotic form is known, the matching procedure can be used to constrain 
the coordinate transformation taking corotating coordinates {t,x,y,z) to internal coordinates {T, X,Y, Z). This is 
the approach to matching taken in this paper. In the next two sections, I calculate the second black hole's tidal fields 
and the perturbations they induce on the first black hole. 

A. Tidal fields of the companion black hole 

Thorne and Hartle |^ have analyzed the motion of an isolated black hole in an arbitrary surrounding spacetime. 
They define and discuss the black hole's local asymptotic rest frame (LARF). In the LARF of BHl, the metric can 
be expanded in powers of the black hole's mass toi as follows^: 

g = gt"' + mi,g(i) + TO2g(2) _|. . . . _ (3 3) 

Here the metric g^'^^ represents the external universe without BHl; the rest of the terms represent the black hole's 
internal gravitational field and the nonlinear interaction between internal and external fields. In this section, I will 



^This expansion is written in Eq. (2.5) of Q; I have substituted ttii for M in that equation. 
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focus on the external metric g^^^ and use it to constrain internal perturbations. Throughout this section and in the 
rest of this paper, boldface letters denote spacetime tensors of all ranks (including 4- vectors). 

In the case of BBHs, the external metric is simply that of a single black hole; it is the metric of the companion 
black hole BH2 of mass m2- This metric must be expressed in LARF coordinates. With this goal in mind, consider 
first a freely- falling observer in a circular, equatorial orbit around a Kerr black hole of mass 1712- (I will later specialize 
to a non-rotating black hole.) The Kerr black hole represents BH2 while the observer's local Lorentz frame (same as 
proper reference frame) represents the LARF of BHl. 

The metric near the observer's world line is determined by the Kerr black hole's electric-type and magnetic-type 
tidal fields as seen in the observer's local Lorentz frame. These tidal fields can be evaluated by taking components of 
the Kerr spacetime's Weyl tensor C in a parallel-propagated orthonormal (PPON) tetrad along the observer's world 
line. The vectors in this tetrad form the coordinate basis of the local Lorentz frame at the location of the geodesic 
orbit. 

The electric- type tidal field as seen by such an observer has been calculated by Fishbone |^ and Marck . Marck 
has computed a PPON tetrad (Aq, Ai, A2, A3) along arbitrary geodesies of the Kerr spacetime, with Aq equal to the 
4-velocity of the geodesic. He obtains the electric-type tidal field by evaluating 



R 



OiOj 



C(Ao, Ai, Aq, Xj 



(3.4) 



I specialize his tetrad to circular, equatorial geodesies; I also label the tetrad vectors (and hence coordinate axes) 
differently. Initially (that is, at proper time T = 0), I choose Ai to be radially outward (in Boyer-Lindquist coordi- 
nates); A2 is chosen so that the projections of Aq and A2 on a constant-Boyer-Lindquist-time-i surface are parallel; 
and A3 is then chosen to give (Ai, A2, A3) positive (i.e., right-handed) orientation. With this choice of tetrad, I obtain 
the magnetic-type tidal field using Marck's work by evaluating 



R 



Oijk 



C(Ao, Ai, Aj, Xk 



(3.5) 



The results of the calculations (3.4) and (3.5) with the above choice of tetrad are 



i?oioi = ^ ^1 - 3 (^1 + ^ j cos' 

-K0303 - I + ) ' -^^0102 - 

^0112 = ^^0121 = -R0323 = — -R0332 ~ 

-R0212 — —R0221 = Rma.1 ~ ^^0313 ~ 



- Ro2Ql 

3to2>V 
3to2>V 

d4 



^0202 — 

_ 3m2 
n;2 



1 



1-31 



sin' rjT 



cos riT sin riT, 



1/2 



cos ilT, 



1 



1/2 



sin fiT, 



(3.6) 



and the rest of the Weyl-tensor components are zero. Here T is proper time along the geodesic, 7712 is the mass of the 
Kerr black hole, d is the Boyer-Lindquist radial coordinate of the circular, equatorial orbit, and 

(3.7) 



is the (exact) rotation rate of the black hole's tidal field as seen in the local Lorentz frame. The quantity W is given 
by 



d 



\/m2/d± a/d 



1 — 3m2/d =F 2ay^m2/d^ 



1/2 ' 



(3.8) 



where maa is the black hole's angular momentum. The upper sign in Eq. ( |3.8| ) is for a retrograde orbit while the 
lower one is for a direct or prograde orbit. 

Notice that for d 3> m2, the electric- and magnetic- type tidal field components at T = are simply related via a 
Lorentz boost with low velocity (1712/ dy^^. For example. 



i?0112|^^0= -(m2/d)^/2[i?0101 +i?1212]^^0 = -(m2/d)l/2[2i?0101 +^0202]y^0 

to lowest order in 7712/0?. This fact will be used later in this section. 



(3.9) 
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In the local Lorentz frame, the spacetime metric can be written as a n exp ansi on in powers of distance TZ from 
the observer's geodesic world line [Q. The two types of tidal field {\A) and (3.5) determine the metric up to and 
including terms of 0{TZ'^). After exploiting some gauge freedom (see Sec. V.A.2 of [|12| and Eq. (2.7) of [||), the metric 
can be written in local coordinates {T,X,y,Z) as 

500 = -1 - Ra^03{r)X'X' + 0(11^), 



9oi 



S^J [1 - i?0fe0m(r)A'"A""] + 0(7^3), 



(3.10) 



where TZ= {X^ + 2^)1/2 Substituting the expressions ( |3.6D in Eq. ( |3.10| ) gives 

, ma L . 3W2 

500 = -1 + ^ 



9ox 
9oy 
9oz 
9^j 



2m2>V 
2m2W 



(P 



{{Z^ - 3^2) sin m - Xy cos vlt] , 



1 + — j (3^cosi7r- A'sinr2r)z, 



3(1 + ^) {X cos^T + ysmnrf 



dp 



(3.11) 



up to and including terms of 0{TZ^). 

The rotation rate is only correct for test-particle orbits and is exact in that case. The correct rotation rate 
of the second black hole's tidal field — measured in a local inertial frame in the first black hole's LARF — is actually 



determined b y the post-Newtonian metric (2.9) and by the requirements that (i) this metric match the LARF metric 
[given in Eq. ( 3.24 ) in Sec. Ill B below]; and (ii) th e LAR F coordinate system be non-rotating relative to local inertial 
frames. The rotation rate is calculated in Sec. [VA by transforming the metric ( p.9[ ) to internal coordinates and 
requiring a match to the LARF metric ( 3.24 ). There it will be seen that the rotation rate is^ 



1-^+0(6^) 



(3.12) 



Note that post-Newtonian corrections to the orbital angular velocity u! of 0{e^uj) have not been included in th is pap er 



to 



The metric ( 3.11 ) is valid for all radii d which allow a circular, equatorial, geodesic orbit. To apply Eq. ( 3.11 
the situation of widely-separated non-rotating BBHs, I specialize to a Schwarzschild black hole by setting a = and 
take the limit of small ma/d, keeping only lowest-order terms in 1112/ d. [In particular, I replace W/d with (ma/d)^/^.] 
I then replace d with b, C l with il, and local coordinates (T, X ,y,Z) with internal coordinates (T, X, F, Z) [which are 
described above Eq. (3J)]. The result is 



ffoo = -1 + ^ [3(Xcosr2r + ysinOr) 



2 R^] 



9ox 

goY 
9oz 



[{Z^ - Y^) sin - XY cos flT] , 
[{X^ - Z^) COS riT + XY sin QT] , 
(Fcos riT-X sin flT)Z, 



2m2 /to2 
2m2 /m2 
2m2 frn2 



|l + ^ [3{X cosnr + Y sinnrf - R'^]^ 



(3.13) 



This rotation rate can also be calculated by looking at geodetic precession of parallel-propagated vectors in the LARF ll3 
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where R= {X^ + Z^)^/^ as defined in Eq. (U). 

The metric ([3.13| ) is still not applicable to BBHs since the observer was taken to be massless [there are no f actors 
of mi in Eq. ( 3.13| )]. This can be fixed easily. As mentioned above, the factors of (7712/6) in goi in Eq. ( 3.13| ) arise 
from a Lorentz boost with low velocity (7712/6)^/^. But the correct (Newtonian) relative velocity between the black 
holes is e = [(7711 + 7772 )/6]^^^. So I replace the factors of (7772/6)^/^ in Eq. ( |3.13 ) by (777/6)^/^. The resulting metric 
includes the second black hole's tidal fields but does not include the first black hole's gravitational field: 



500 = -1 + ^ [3{X cos ilT + Y sin QT) 



2 R^] 



gox 
goY 
90Z 



2to2 pm 



27772 prn 



J — [{Z^ - Y^) sin nr - XY COS OT] 
[{X^ - Z^) cos nr + XY sin nT] 



2 777 o / 771 

— ^ J —(Ycosnr ^ X sin nr) z, 

0-^ V 



= 5,^. 1 1 + !^ [s{x COS nr + Y sin nr)^ - } • 

In the buffer zone around BHl, this metric provides the asymptotic form of the perturbation on BHl. 



(3.14) 



B. Schwarzschild perturbation 



The next stage is to solve the LEEs abo ut th e Schwarzschild metric for a perturbation which is finite at the horizon 
R = 7711/2 and asymptotes to the form ( |3.14[ ) as R/mi 00. For ease in dealing with the LEEs, I transform to 
spherical, isotropic, internal coordinates (T, R, 6, (p) by letting 



X — R sin 9 cos < 



Y = R sin 9 sin ( 



The unperturbed Schwarzschild metric in these coordinates is 

\+'^)\dR^ + R\d9^ 



Rcos9. 



mi/2R" 



- 777l/2i?^ 



dT^ 



- sin^ 9 d(j)^)] . 



(3.15) 



(3.16) 



The metric (3.14) in these coordinates is 



'dT^ + dR^ + R^{d9'^ + s\T?9d(j)'^) 

7712-R 



„ .i—dT [cos 9 sin(0 - Q.T)d9 + sin 9 cosl29) cosU - Q.T)d(j3\ 
b-^ V 



m2R? 



[3sin2 6'cos2(0 - VlT) - l] [dT^ + dR^ + R^{d9'^ + sin^ 6* d^^)] 



(3.17) 



The linearity of the LEEs allows me to look separately for solutions corresponding to the electric-type and magnetic- 
type tidal fields. First I look for a perturbation hi = g — gg (where gs is the unperturbed Schwarzschild metric and 
g is the full metric including the perturbation) of BHl which corresponds to the electric-type tidal field of BH2 and 
is of the form 



hi 



7772 i?^ 



[3sin2 6icos2(0 - VlT) - l] [h{R)dT'^ + f2{R)dR^ + h{R)R^{d9'^ + sin^ 6I d./-^)] , 



(3.18) 



as suggested by Eq. ( |3.17 ). In this notation, dT, c?i?, d9, and dcj) are coordinate one-forms and dT^ denotes the tensor 
product dT®dT. The functions /i, /2, and /a are to be determined by solving the LEEs with the following boundary 
conditions: (i) /i(i?), f2{R)i and faiR) ar e req uired to approach 1 as R/nii — > 00 so that the perturbation ( 3.1§| ) 
matches the electric-type tidal field in Eq. ( 3.17 ); and (ii) hi is required to be finite at i? = 1711/2. 

Consider solving the LEEs order by order in e = {m/bY^^. Time derivatives of the components of hi produce factors 
of 777.1 ~ e"^ in the LEEs and can thus be neglected. A solution for hi can then be found using the Regge- Wheeler 



formalism |14| for analysis of stationary Schwarzschild perturbations. Regge and Wheeler decompose perturbations 
into even- and odd-parity modes and analyze them in a particular gauge chosen to simplify computations. In their 
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classification hi is a superposition of static^ even-parity modes with angular numbers I = 2 and m = —2, 0, 2. The 
general solution of the LEEs for static even-parity modes with Z > 2 is well-known in Schwarzschild coordinates and 
is given in Sec. IV of ||l5|, for example. A particular solution with I — 2 that is finite at the black hole's horizon and 
contains an arbitrary multiplicative constant is easily obtained from the general solution, and is given in Eqs. (6.5) 
and (6.7) of for example^ After transforming this solution to isotropic coordinates, the multiplicative constant 
is determined by imposing the boundary condition (i) (given at the end of the previous paragraph). This yields the 
following solution for the radial factors fi{R), f2{R), and fsiR) in isotropic coordinates: 



/i(i?) = (l 
f2{R) = (i 

MR) = fi 



2RJ 



2r) \ 
2r) 



1 



mi \ ° 
2r) ' 

2rJ 



2m{ 



(3.19) 



Next I look for a perturbation h2 = g — g,, 
the form 



of BHl corresponding to the magnetic-type tidal field of BH2 and of 



4m2 



R^F{R)dT [cos e sin{(j) - nT)de + sin 9 cos 26 cos(0 - Q.T)d(j)] 



(3.20) 



as suggested by Eq. ( 3.17 ). The function F is to be determined by solvin g the LEEs with the following boundary 
conditions: (i) F{R) — > 1 as R/nii ^ oo so that the perturbation (3.20) matches the magnetic-type tidal field in 
Eq. ( |3.17| ); and (ii) h2 finite at i? = mi/2. As was done for hi, time dependence is ignored in h2 since time derivatives 
produce higher-order terms. In the Regge- Wheeler classification, h2 is a superposition of stationary odd-parity modes 
with angular numbers I ~ 2 and m = —1,1. The general solution of the LEEs for stationary odd-parity modes that 
are finite at the horizon and have Z > 2 is given in Schwarzschild coordinates in Eq. (38) of [Q. This solution is only 
determined up to a multiplicative constant. The particular case Z = 2 is easily obtained from the general solution, and 
is given in Eq. (6.10) of g, for example^. After transforming this solution to isotro pic co ordinates, the multiplicative 
constant is determined by imposing the boundary condition (i) [given below Eq. (3.20)]. This yields the following 
solution for the radial factor F{R) in isotropic coordinates: 



FiR) 



miV 
2RJ 



mi \ ' 
2r) 



(3.21) 



The metric i n the inter nal re gion near BHl is now complete. It is giv en by the Schwarzschild metric (3.16) plus the 
perturbations (3.18) and ( 3.20| ) with radial factors given in Eqs. (3.19) and (3.21); in other words, g = gs + hi -I- h2. 
In spherical isotropic coordinates (T, R, 9, 0), this internal metric is 



ds^ 



R^dT [cos 9 sin(0 - VLT)d9 + sin 9 cos(26l) cos(0 - nT)d4)\ 



4m 



\l b V 2RJ V ' 2RJ 
m2R^ 



63 



[3 sin^ 9 cos^ i<^-m-l]\(l-^)dT^'+{l-^) {l + ^)dR^ 



miV 
2r) 



(1+2]^) 



2ml 



i?2 



R^{de^ +sm^ 9 ( 



(3.22) 



In isotropic coordinates (T, X, Y, Z) , this metric is 



'^Time dependence in Eq. (3.18) is to be ignored, as explained above. 

^The notation in Sec. VI of may be confusing: R there denotes a dimensionless quantity obtained from the Schwarzschild 
radial coordinate ra hy R — Vs/M where M — in my notation mi — is the mass of the black hole being perturbed. 
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500 = 
9ox 

90Y 

9oz 
9ij 



^{'-^fi'^^y [(^^ - ^^^^ - ^^-^^^] ' 

(l - ' (l + ' [iX' Z^) COS + XY sin QT] , 
{Y cosnr- X sin nT)Z, 



2r7i2 [rn 
2m2 / »n 



63 V b 2i?J ' 2RJ 



l + ^y + ^ cos l^T + y sin r!T)2 - R^] 



mi\4 2to? 

IH -] 

2R) i?2 



2mi I iii,-^ 
On ^ 1 



R 



TO? \ X'X^ 



4i?2 ; i?2 



(3.23) 



Expanding the components (3.23) in positive powers of rrii/R and R/b in the buffer zone nii <^ R b and keeping 
only lowest-order terms yields the LARF metric: 



500 
9ox 

90Y 

90Z 



2toi Tot r , „ „ on 

-1 + — i + -i- \'3(x cosnr + Y sinnry - r^] , 

R 0-^ 

'^^[(Z^^Y^)smnT-XYcosnT] , 
^ Ux^ - Z^)cosnT + XY sin QT] , 



2m2 /to 



(y cos i7r - X sin nr) z, 



9ij = S^j i 1 + 



2toi TO2 



[3{x cos r^T + r sin nT)^ - I 



R 63 

This metric includes the first black hole's (weak) gravitational field as well as the second black hole's tidal fields 



(3.24) 



IV. DISTORTED-BLACK-HOLE METRICS IN COROTATING COORDINATES 

The post-Newtoni an rn etric (^.91), when expressed in internal coordinates (T, X, Y, Z) in the buffer zone around BHl, 
must take the form (3.24). The next step is to find explicitly the coordinate transformation in the buffer zone taking 
corotating post-New tonia n coordinates to these internal coordinates. Applying the inverse of this transformation to 
the internal metric ( |3.23| ) will put that metric in corotating coordinates (t,x,y,z). An identical procedure will then 
be followed to obtain the metric near BH2 in corotating coordinates. 



A. Buffer-zone coordinate transformation 



In this section, a serie s of c oordinate transformations are performed on the metric (2.9) in the buffer zone of BHl 
to bring it to the form ( 3.24 ). Composing these transformations gives the final transformation from corotating to 
internal coordinates. Throughout this process terms of 0{m?) are dropped; justification for this will be given at the 
end of the section. 

Begin with the near-zone metric (2.9) with terms of 0{m?') removed. Restrict attention to the buffer zone TOi <C 



ri <C 6 since this is where the corotating coordinate system and internal coordinate system overlap. Center the 
coordinate grid on BHl by shifting the origin to (x, y, z) = (to26/to, 0, 0). This is done by defining a new coordinate 



TO26 



(4.1) 



Next expand the metric in powers of the distance 
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ri = (^2^y2^z2)l/2 (4 2) 

to the new origin. The expansion for is 

-=[{b + 0+y+z] =i-^ + ^3 + ■••• (4-3) 

Positive powers of ri in the metric components come in the form {ri/b)P with integer p > 0. Since ri ^ 6 in the 
buffer zone, discard terms of 0[{ri/b)^] or higher. This results in the foUowing metric: 

ds^ ^ dt^ [-1 + 2mi/ri + (2to2/6)(1 + m2/2m) + (m2/fe^)(2£2 - - z^) + {mlb^){(^ + y^)] 

~ 2ujydtd£, [l + 2m2/5 + 2?7ii/ri - 2m2£/6^] + 2Ludtdy^^{m2/'m){b + Ami + 2m2) - 2fib/ri 

+ C(l + 2mi/ri - 2n/b) + im2/b^) [2mi^^/m - (1 + mi/m){y^ + z^)] | 

+ {d^^ + dy^ + dz^) [1 + 2mi/ri + 2m2/b - 2m2^/9 + {m2/b^){2^'^ - y^ - z^)] . (4.4) 
Now renormahze the time-coordinate by defining 

i = t[l + (m2/6)(l + m2/2m)], (4.5) 
and then perform a partial Lorentz transformation by setting 

t = i + {m2Uj/m){b + 4mi + 3m2 + rn^/2m)y, 

£,=x, y = y, z = z. (4.6) 



In the new coordinates {t,x,y, z), the metric (4^) is 



ds^ = di [-1 + 2mi/f + (m2/6^)(3i2 - f^) + {m/b^){S:'^ + f)] 

— 2ujdidx {y [1 + (m2/2m6)(6r7ii + 7to2) + 2mi/f] — 2m2xy/b'^} 

+ 2ujdidy {x [1 + (m2/2m6)(-2mi + 3m2) + 2mi/f] + (m2/6^)(3x^ '2z^)} 

+ [di^ + df + dz^) [l + 2mi/f + 2m2/5- 2m25/6^ + (m2/fe^)(3i^ - f^)] 

+ {m2/b)dy [{m2/m + 2i/b)dy - {2y/b)dx\ (4.7) 

where f = (i^ + y^ + z^Y^^ and terms of 0{m?) have been dropped, as is done throughout this section. 
Next clean up the spatial part of the metric by putting 

i^i, x^ x{l - m2/b) + (to2/26^)(£^ +f- z^), 

y = y[l — (m2/2TO6)(2mi + 3m2)], z — z{l — m2/b) + {m2/b'^)xz. (4-8) 



Transforming the metric (4.7) using Eq. ( [4.q ) results in 

ds^ = dt [-1 + 2mi/f + (m2/&^)(3x^ - f^) + {m/b^){x^ + f)] - 2uodidx [y(l + TO2/& + 2mi/f) - TO2iy/6^] 
+ 2udidy [x [1 ~ (m2/m6)(3mi + 7712) + 2mi/f] + (m2/2fe^)(7i^ - 3?/^ - 55^)} + 2{m2Uj/b'^)yzdidz 
+ (di^ + df + dz^) [l + 2mi/f + (to2/&^)(3£^ - r^)] (4.9) 

where f = {x^ + f + z^Yl"^ . 

Focus attention on the terms 2u!di{xdy [I — {m2/mb){3mi + 1712) + 2mi/f] — ydx{\ + TO2/6 + 2mi/f)} in Eq. (ff.S|). 
These terms contain information about the rotation of the coordinate axes. However, they are not yet in the form 
of the rotation terms 2Q.{1 + 2mi/f)dt{xdy — ydx) that result from rotating — at a constant rate Vl and in an active 

„2 

sense, i.e., using a pull-back map — the metric ds^ = dt (—1 -I- 2mi/f) -|- [dx^ + df + dz^){l + 2mi/f), which is a 
fragment of Eq. (4.9). An additional coordinate transformation is required to bring the former terms into the latter 
form. With this goal in mind, look first for a gauge transformation taking the perturbation 

7 = 2ujdt{xdy [1 — {'m2/mb){3m.i + TO2) + 2TOi/f] — ydx{l + m2/b + 2TOi/f)} (4-10) 

on a flat background metric ds^ = —dP + dx^ + df + dz? to the perturbation 
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7 = 2f7(l + 2mi/f)dt{xdy — ydx). 
In other words, look for a vector field rj such that 



(4.11) 



(4.12) 



In order to solve Eq. (1.12 



perturbations ( 4.1C ) and ( 4.11 



it suffices to co nsider J7 with only one nonzero component 77* = rf{x,y,z). The 
) when put in Eq. ( 4.12| ) yield 



- ^y{l + 2mi/f) = = 7^^ + dif /dx — —ujy{l + 7712/6 + 2rni/f) + dif /dx, 
Q,x{l + 2?Tii/f) = 7j-g = 7j-g + drf I'dy ~ lux [1 — {m.2 / mb){im.i + 7712) + 2?7ii/f] + drf jdy. 



(4.13) 
(4.14) 



These two equations determine the rotation rate as follows: the function 7/*(x,y, z) must satisfy d^rf / dxdy = 
d^rf jdydx. Taking djdy of Eq. ( 4.15 ) and djdx of Eq. ( 1.14 ), equating the mixed partials of 77*, and ignoring terms 
of 0(m}) yields the following equation for fi: 



f2 - w [1 - (m2/6)(l + 2mi/TO)] = uj(\ + 771,2/6) - 17, 



(4.15) 



which has solution = a; [l — /i/6 + 0(e"^)] . This is the rotation rate o f the sec ond b lack hole's tidal field as s een in 
the fi rst b lack hole's LARF; this value confirms the claim in Sec. Ill A [see Eq. ^A%\. With f7 in hand, Eqs. (|4T|) 
and (4.14) now yield 



jf — {m20j/b){l + mi/m)xy. 



(4.16) 



Gauge transformations can also be thought of as resulting from infinitesimal coordinate transformations. The 
coordinate transformation corresponding to the gauge transformation given in Eqs. (4.12) and (4.16) is 



t — t — {m2L0 /b){l + Tni/m)xy, 
x^x, y^y, z = z. 

The metric (|4.9|) expressed in the new coordinates (t, x, y, z) is 



(4.17) 



ds^ = df [~1 + 2mi/f + {m2/b^){3x^ - f^) + {m/b^){x^ + f )] + 2f7(l + 2mi/f)di{xdy - ydx) 
+ {m2Uj/b'^)dt [2xydx + (7x^ - - 5z^)dy + 2yzdz] 
+ {dx^ + df + dz^) [1 + 2TOi/r + (7772/6^) (3^2 _ f^)^ ^ 



(4.18) 



where r = (x^ + + z^y^^. 

The next step is to undo the rotation of the coordinate system. But first some fine-adjustment of coordinates is 
needed in order to obtain the LARF metric (3.24). To find out what is required, the metric (3.24) can be put in 
coordinates rotating with angular velocity fi. It turns out that the fine-adjustment needed is 



i = t + {m2Uj/2b^){3x^ -y - z^)y, 
x = x, y = y, z^z. 



(4.19) 



In the new coordinates {t,x,y, z), the metric (1.18) is 



ds^ = dr 



+ 2mi/f + (7r72/6^)(3x2 - f^) -|- {m/b^){x^ + f)] + 2f](l + 2mi/f)dt{xdy - ydx) 



+ 4{m2Lu/b'^)dt [-xydx + {x? - z'^)dy + yzdz] + {dx"^ + df + dz"^) [l -I- 2mi/f + (m2/6^)(3x^ - r^)] , (4.20) 
where f = {x^ + f + z^)^/^. Now eliminate the rotation of coordinates by defining 



t = T, X ^ X cosnr + Y sin nT, 
y ^ -Xsinnr + YcosnT, z ^ Z. 



(4.21) 



Transforming the metric (1.2C) using Eq. (1.21) results in the LARF metric (3.24). 

The transformation from corotating post-Newtonian coordinates {t,x,y,z) to isotropic internal coordinates 
(r,X,r, Z) can now be gotten by composing the transformations (|t]), (U), (^), (Q, ( [1.17D , ( |4.19| ), and ( [1.21D . 
Inverting this composite map gives the following transformation from internal to corotating coordinates: 



15 



T = 




X 
Y 

Z 



rcos nr - AsinQT, 
rsinf7T + Acos ftT, 

~b W 



1 



(4.22) 



where 



A = y 



m 
T 



m2 
2m 



b J 



= uj[l- 



and terms of 0{m^) have been dropped. In terms of the coordinates {x,y,z), 



i?-(r2 



Z2)l/2 



(4.23) 



(4.24) 



[cf. Eq. ^)]. 

There are two reasons why terms of 0{m^) were dropped from the metric ( |2.9| ) at the beginning of this section. 
First, suppose that such terms were kept and were used to calculate higher-order deformation of the black hole. Since 
internal metric components are coupled to each other via the Einstein equations (in particular, the components of a 
black hole perturbation are coupled via the LEEs), to be fully consistent, terms of O(to^) would have to be included 
in the spatial part gij of the metric (2.9). But these terms are of higher order than first post-Newtonian, and so have 
not been included in this paper. 

Second, black-hole perturbations with asymptotic form m'^r^/b^^'^ {p > 2) in the buffer zone, which come from 



terms of 0{m'^) in qqq in Eq. (2.£), are actually smaller in the internal region than the perturbation with asymptotic 
form m2ri/b'^ in the buffer zone; the latter perturbation has been ignored in this paper. Once terms of Oim?) were 
dropped in Eq. (^^), all terms of 0{m^) were consistently discarded in this section. 



B. Internal metric in corotating coordinates 



In this section, 
in the buffer zone 
preserve finiteness 
when performing 
first defining = 
transformation (4 
coordinates (T, F, 
functions Pap are 



the transformation (4.22) is applied to the internal metric (|3.23 ) throughout region I (not just 
). This puts the inter nal m e tric i n cor otatin g post- Newtonian coordinates (t,x,y,z). In order to 
of the perturbations (3.18)-( ^.1£ ) and ( 3.20|)-( 3.21) at the horizon of BHl, all terms must be kept 
the transformation. The rotation in Eq. ( fl.22 ) can easily be performed on the metric (3.22) by 
(p — riT and then setting F = i? sin 6* cos A = RsinO sinip, and Z = Rcos 9. To complete the 
22), define the functions Pa/3{x,y, z) for a,P — 0, ..,3 to be components of the i nternal metric in 
A,Z)\ write the components as functions of {x,y,z) using Eqs. (4.22) and ( 4.23| ). Explicitly, the 



Pr 



'00 



—, — +^(1 -] 3F2-i?2 ^f7F(l -] {1 + — ] (2Z^ ~ R^) 

l + mi/2i?y 63 V 2rJ ^ ' b^ \ 2rJ V 2rJ ^ ' 

2ml'' 



Poi = Pi 



P02 = P 



20 



2RJ (r^ + A^) i + ^(3r^ 



r' 



V 2rJ 



i?2 



2em2 

P~ 2RJ ' 2RJ 
2em2 miNj2 (^^m^- ^2 



1-^) 
2rJ 



(3F2 - R^] 

2 / -m, \ 4 



/ mi\4 2m: 



2RJ 



2m^ 



i?2 



16 



P03 - 



13 



30 



Pll = 1 



P22 = 1 



33 



2em2 

2i?y 

mi \ ' 
2r) 



1 



P\2 — P2I — — 



2RJ 

2r7iim2 



P31 



P23 — P?,2 — 



2r7iim2 
2r7iim2 

i?353 



mi \ ^ 
" 2^/ 



2i?y 



\4 2mf 



miy 
2r) 



miy 
2r) 
miy 
2i?/ 



/ mi 



V 2i?y 



i?2 

2mf 
~W 
2ml 
~W 

R^)TZ, 
R^)AZ, 





/ 


2 


i?3 




2mi 


(- 


mj 


i?3 




2mi 


(- 


9 

mj 


i?3 





(4.25) 



where e = (m/6)^/^, = a;(l — -R — (r^ + A^ + Z^)^/^, and F, A, and Z are given in terms of {x,y,z) in 

Eqs. and p^ ). 

Next define the functions K^{x,y,z) for p, ct = 0,..,3 by is:° = dT/dxP, = 9r/aa;'', = dA/dxP, and 
if3 = dZ/dxP, where (T, r,A, Z) are to be expressed in terms of (t,x,y,z) using Eqs. (4.22) and ( 4.23| ). Exphcitly, 



the functions are 
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m2 
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1 _ m2Z 
i^? - 0, 



m2 _ m^S, 
~b 6^' 



(4.26) 



where e = {m/by/^ and ^ = a: — m2b/m. 

The metric in region I can now be written in terms of the functions Pq,^ and Kp . It is given in corotating coordinates 
{t,x,y,z) by 



{x,y,z)= ^ Paa{x,y,z)K'^{x,y,z)K^{x,y,z) 



(4.27) 



with Paf3 and if^ as defined in Eqs. (4.25) and (4.26). Note that the metric components are exphcitly independent of 
time t. This metric is vahd throughout region I (up to the first black hole's apparent horizon) and matches (to some 
finite order; see below) the post -New tonian metric (|]^) at ri = r™. 

Errors in the internal metric ( 4.27 ) will only be analyzed in the wea k-gravity buffer zone mi <^ ri <^ b. The largest 
errors come from inaccuracies in the coordinate transformation (4.22). Terms of the form [m? /b'^){ri/bY for integer 
p> 1 have been ignored in Eq. ( 4.22 ). This leads to errors Sg/^iy ^ for ri ^ r™. 

The match between the internal and post-Newtonian metrics at ri = r™ is not perfect; there are discontinuities 
[dfiiA ill the metric components on that 3-surface. A term of the form mf/rf in the internal metric component goo [^^ 
given in Eq. (4.27)] is not matched in the post-Newtonian metric component goo in Eq. (2.9); as a result, [goo 



£3 at 



ri 



Similarly, a term of the form mi/rj is not matched in gij, so [gij] ~ e . Lastly, a term of the form m\/ R 



in the internal-coordinate metric component 500 in Eq- (3.23) gives via a (partial) Lorentz boost an u nmatched term 
of the form emf/rl in the internal-metric component goi in corotating coordinates [given in Eq. ( 4.27| )]; so [goi] ~ e'^. 

The internal metric (4.27) contains terms of the form {Tni/ri)P{m2rf/b^), p > 1, in the buffer zone. These terms 
represent the first black hole's multipole moments and the nonlinear interaction of internal and external gravitational 
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fields. They are of 0(6^+^), p > 1, in the buffer zone and have not been matc hed to the post-Newtonian near-zone 
metric (2^). At the level of accuracy achieved in this paper, the metric (2^) need not be modified to include the 
near-zone gravitational effects of the black holes' deformation. 



C. Metric near the second black hole 



An identical procedure can now be followed to obtain the metric in corotating post-Newtonian coordinates in 
region II. However, it is not nece ssary to repeat all the steps. Thi s me tric can simply be gotten as follow s: exc hange 
TOi <-^- 7712 in the internal metric (3.23) and in the transformation ( [4.22[ ); take x — > —x and y — > —y in Eq. ( 4.22| ); then 
transform. In other words, the metric components in region II [denoted gjji^{t, x, y, z)] are related to those in region I 
[gfi^{t, X, y, z)] by gjj,y{t, x, y, z) — (— l)^g^^(i, — x, — y, z)(with mi ^ ^,2), where p is the number of the indices /i and 
V that are equal to 1 or 2. 

Define Pap to be Pap with mi and m2 exchanged, and similarly to be with mi ^ m2. Then the metric in 
region II is given in corotating coordinates (t, x, y, z) by 



9tiiy{.x,y,z) ^ (-l)P 



3 

E 

a.cr— 



Pa.(T ( 



^V, z)K"{~x, -y, z)K1{-x, -y, z) 



(4.28) 



where p is, as above, the number of the indices and v that are equal to 1 or 2. Again, the metric components are 
explicitly independent of time t. This metr ic is valid up to the second black hole's apparent horizon and matches (to 
a finite order) the post- Newtonian metric (|2.9| ) at r2 = r™. Error analysis for this metric is identical to the analysis 
above for the metric in region I. 



V. RESULTS AND DISCUSSION 



The result of this calculation is an approximate solution to Einstein's equations representing two widely-separated 
non-rotating black holes in a circular orbit. The metric has been expressed in a single set of coordinates valid up to the 
black holes' apparent horizons; the coordinate system chosen is corotating coordinates (i, x, y, z). In th ese co ordinates, 
the metric comp onents are explicitly in dep endent of time t. The metric is sp ecified in region I by Eq. ( 4.27 ), in region 



'2 : 



II by Eq. ( 4.28| ), in region III by Eq. (2.£), and in region IV by Eq. ( ^.1^ ). At the boundaries ri — rj'", r2 
and r = r°"' of these regions, there are discontinuities in the metric components that result from truncation of 
perturbative expansions and finite-order matching. The magnitudes of these discontinuities can be reduced by taking 
this calculation to higher orders. 

The full 4-metric is summarized below: 
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°'2 

^xy + i (e^iV - 625-) sin 2tjr 

+ — (x cos ujr ~ y sin wr) cos 2a;r, 



in region I; 
in region II; 
in region III; 
in region IV. 
in region I; 

in region II; 
in region III; 



Eu,,.=o P^i^^^ y^ z)K'^{x, y, z)K^{x, y, z), 



0, 



0, 



in region IV. 
in region I; 

Tf^,.u=o Pf^i'i-^^ -y, z)K'^{-x, -y, z)K^{-x, -y, z), in region II; 
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in region I; 
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in region IV. 



-xz, 



Y.n,u=o P^^v{x, y, z)K^{x, y, z)K^{x, y, z), 
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-yz, 
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(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 
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In the expressions above, m — mi + TO2, M — ™im,2/m, 5m = mi — m2, e = {m/b) 



1/2 



(m/b^) 



1/2 



[{x — m2b/m)^ + y^ + z^Y^"^, and r2 = [(a: + mib/m)^ + y^ + z^]^/^. Region I is specified by ri < {niiby^^ and region 
II by r2 < {m2b)^/'^ (but these regions do not extend inside the black holes' apparent horizons). Region III is specified 
by Ti > (miby '^, r2 > {m2by/^, and r < b/2e; and region IV by r > b /2e. The fu nctio ns A, B, D, E, N, and S are 
defined in Eq. ( ^.191) . The functions Paf) and K'^ are defined in Eqs. ( [4.25D and (|^). The funct ions Pap and K'^ 
are obtained from P^^ and respectively by exchanging mi and m2. 

The errors and discontinuities in the metric components are summarized in Table |[ The discontinuities should 
be smoothed out before initial data are extracted from the metric. In addition, initial data taken should be relaxed 
numerically to approach more closely an exact solution of the constraint equations. It is expected that higher-order 
versions of this calculation will differ by smaller amounts from an exact solution of the Einstein equations. 
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TABLE I. Errors and discontinuities in the metric components in corotating coordinates. Numbers denote orders in 
e = e.g., 4 denotes 0{e*). The last two columns contain normalized errors. 
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